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Abstract. We study a class of zero-range processes in which the real-space 
condensation phenomenon does not occur and is replaced by a saturated condensation: 
that is, an extensive number of finite-size "condensates" in the steady state. We 
determine the conditions under which this occurs, and investigate the dynamics 
of relaxation to the steady state. We identify two stages: a rapid initial 
growth of condensates followed by a slow process of activated evaporation and 
condensation. Wc analyze these noncquilibrium dynamics with a combination of mean- 
field approximations, first-passage time calculations and a fiuctuation-dissipation type 
approach. 
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1. Introduction 

Real-space condensation, whereby a finite fraction of a system's mass accumulates 
within a microscopic region, is a spectacular phenomenon that is observed in a wide 
range of dynamical systems. For example, it is manifested experimentally in shaken 
granular gases in which particles can diffuse between compartments [1, 2]: as the driving 
strength is reduced, the sand grains cluster into a single compartment. One can also find 
examples in models of macroeconomics, whereby a large fraction of the available wealth 
is accumulated by a single individual [3], and of traffic flow, in which buses serving a 
single route cluster together [4]. 

One of the requirements for a thermodynamic condensation transition to occur as 
the total density of particles is increased is the absence of any restriction on the total 
number of particles that may occupy a single site. In this work, we are interested 
in the case where such unbounded growth of particle number is inhibited. This can 
happen quite naturally within specific applications: for example, the compartments in 
the granular gas experiments of [1, 2] are of finite size, and once they contain more than 
a certain number of particles, any extra particles may diffuse freely out of them. As 
we will show below, a vestige of the condensation phenomenon may still be observed 
in the form of a separation of the system into high- and low-density sites: we call 
this saturated condensation. The key questions then are: (i) Under what conditions 
is saturated condensation observed? (ii) How is the state of saturated condensation 
approached dynamically from some given initial condition? 

In this work, we will provide answers to both these questions with reference to 
a specific well-studied mathematical model, the zero-range process (ZRP, see [5, 6] for 
reviews) . This is a stochastic dynamical model in which particles occupy sites of a lattice 
(or network) and hop to neighbouring sites at a rate that depends only on the number 
of particles on the departure site: it is for this reason that it is described as 'zero-range'. 
The model was introduced and solved for its steady-state behaviour by Spitzer [7], and 
it was realized a little over ten years ago that real-space condensation [8] is possible at 
sufficiently large particle densities under certain conditions on the hop rates [4], even 
within a spatially homogeneous system (see also [5, 6]). The dynamics of condensation 
onto a single site has also been of interest. Here, the focus has been on the late-time 
coarsening of the excess mass into a decreasing number of increasingly massive clusters, 
which in a finite system ends with a process of mass exchange between the last remaining 
clusters [9, 10]. In the steady state, for a finite system, the condensate occasionally melts 
and reforms on a new site: the timescale of this process has been the subject of some 
discussion [10, 11, 12]. 

Variants of the ZRP where mass accumulates on multiple sites — and in particular 
an extensive number of sites — are comparatively little studied. Schwarzkopf et al [13] 
examined the statics and dynamics of a ZRP with transition rates chosen in such a way 
that a single condensate is destabilized in favour of either a finite number of extensive 
condensates, or a subextensive number of subextensive 'mesocondensates'. In both 
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cases, one still has a finite fraction of the mass occupying a vanishingly small fraction of 
the sites in the thermodynamic limit, and thus a true condensation transition is observed 
in this model. 

By contrast, we consider here the case where the hop rates in the ZRP are chosen 
such that the mass of condensates reaches a finite size that does not increase with 
system size. This prevents a true condensation transition, but nevertheless admits 
the possibility of saturated condensation discussed above. Since the true condensation 
limit can be approached by taking the upper limit on the size of a condensate to 
infinity, it seems clear that saturated condensation can have interesting consequences 
at a phenomenological level even if no formal singularity remains. This is similar to 
the equilibrium phenomenon of "micellization" in which attractive particles can form 
clusters whose size is limited by their packing geometry [14] . An analysis of this problem 
as an instance of saturating condensation in equilibrium was offered by Goldstein [15]. 

After recalling the definition of the ZRP in section 2 and briefly reviewing the 
conditions for a condensation transition in the homogeneous system (section 2.1), we 
present in section 2.2 the conditions on the ZRP hop rates for saturated condensation 
to occur. We then turn to the dynamics of the process. One can first ask about 
the dynamics within the steady state. This involves evaporation and formation of 
condensates, and in section 2.3 we calculate the rates at which both processes take 
place. The remainder of our work, presented in section 3, concerns an investigation 
of the relaxation to the steady state from a prescribed initial condition. Our main 
finding is that this is a nontrivial, two-stage process. First, mass rapidly accumulates 
on sites that are selected by the local dynamics and in a way that depends on the initial 
condition. The number of such sites typically differs from its global equilibrium value. 
This gives rise to a much slower second stage in which condensates are nucleated and 
evaporate as activated processes. We obtain a detailed portrait of both stages of the 
relaxational dynamics through a combination of mean-field theory, calculation of first- 
passage properties and an approach reminiscent of a fluctuation-dissipation analysis. In 
certain cases, these approximations agree remarkably well with stochastic simulations. 
These analytical results thus constitute a more complete account of a nonequilibrium 
condensation dynamics that has been achieved so far. 

2. Presentation of the model: Steady state and condensation 

The model we consider is defined on a one- dimensional lattice of L sites with periodic 
boundary conditions. Each site i can be occupied by an arbitrary number of particles 
Ui. Since the system is not connected to any reservoirs, the total number of particles in 
the system = X^j^i constant. A particle can move from site i to a neighbouring 
site i ±1 with rates uf{ni), respectively (sec figure 1). We call vf{n) the hopping rate 
per particle so that uf{ni) — nivf{ni). Qualitatively, our main results hold for both 
symmetric and asymmetric hoping rates but to avoid redundancies we shall only present 
the symmetric case. The definition of the zero-range process (ZRP) is that the hopping 
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rates depend only on the number of particles at the departure site and not, for instance, 
on the occupancy of the target site. In general the rates could also vary from site to 
site but in this work we consider only spatially homogeneous systems. 

mJ(3) «J(3) 



Figure 1. Presentation of the model. The arrows indicate the allowed transitions and 
the rates at which they occur. 



2.1. Factorization and condensation 

Although the ZRP has been extensively reviewed in the literature [5, 6] we shall 
briefly summarize some known results that will be important in our understanding of 
saturated condensation, deferring to these articles for further details. First, because the 
interactions between particles are limited to a single site, the steady-state distribution 
of occupancies factorizes. That is, in the canonical ensemble one has [6] 

P(M) = ^ n 9in^) S ( n, - n) (1) 

' 1=1 j=\ 

where the partition function Z^^^ is given by 

L L 

^^.^ = E n ^(^^) ^ ( E ^» - ^) (2) 

and the factors g{ni) are determined by the hopping rates 

n ^ 

g{n) = JJ — for n > and ^(0) = 1 . (3) 
j=i 

The delta functions in (1) and (2) simply enforce the constraint that the total number 
of particles on the lattice is fixed. The marginal probability that a given site has n 
particles is given by 

Pi{n) = • (4) 

As previously mentioned, the ZRP admits an interesting condensation transition. 
Although it can be worked out directly from the canonical ensemble [16], the condition 
for condensation is most easily seen in the grand canonical ensemble. Introducing a 
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chemical potential /i and the single-site partition function Zi {/i) , the partition function 
for the L-site system reads 

L 



A'"-E"=ilnKi)] 



= [Ml^t (5) 



N 

Note that the only correlations between different sites in the canonical ensemble come 
from the constraint on the total number of particles in (1). Since this constraint has been 
removed in favour of a chemical potential, the sites are now completely uncorrelated in 
the steady state and the L-site partition function reduces to a single-site problem. To 
get more insight into Zi, one can rewrite the hopping rate per particle as v{j) = Vq 
where Vq is the hopping rate of a single particle and e'^'^-'^ encodes the interaction between 
the particles. For instance, h{j) = for all j corresponds to non-interacting particles. 
The partition function then reads 

Zi(/x) = ^exp[-F(n,/.)] (6) 

n 

where F{n, /i) — f{n) — /in and we have introduced 

n 

f{n) = ln{vo) n + ln(n!) + h{j) . (7) 

The marginal probability that a single site has n particles is then given by 

p(n|/x)=e'^"^(n)^=^ (8) 

explpn - f{n)] (9) 



and the average number of particles per site by [6] 

We thus see from (9) that /(n) plays the role of a single-site free energy, and hence that 
F(n, /i) is a single-site grand canonical potential. 

The general idea of ensemble equivalence is to ask what chemical potential fi should 
be imposed to get a given value of (n). To detect a possible condensation transition, 
one thus looks for the maximum density pc for which equation (10) has a solution. If 
this maximum is infinite, then (10) is always solvable and there is no transition. If, 
however, pc is finite then equation (10) cannot be solved for a density greater than pc 
and the excess mass condenses on a single site, which can thus carry a finite fraction of 
the total mass of the system. The breaking of ensemble equivalence between canonical 
and grand canonical ensemble is a signature of the condensation transition. 

The condition for (10) to have a solution has been worked out and yield a criterion 
on the form of u{n) to observe condensation. It can be summarized as follows [6]: 

• if u{n) decays to a non-zero constant more slowly than u{n) ~ /3{1 + 2/n), one 
observes above a non-zero critical density the appearance of a single condensate in 
a background fluid which remains at the critical density. 
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• a u{n) — )■ as n — )■ oo, condensation occurs at all densities and the fraction of 
particles in the fluid phase tends to zero. 

• otherwise, and in particular if u{n) increases as n — > oo, condensation does not 



The first two cases, in which there exists a true thermodynamic phase transition, 
have previously received much attention in the literature (as reviewed in [6]). What 
we shall show in the following is that the third case also may also exhibit interesting 
condensation-like features, despite the absence of a true condensation transition, when 
the stationary state supports a coexistence of sites at two characteristic densities. We 
shall refer to this case as saturated condensation. 

2.2. Criteria for saturated condensation 

As mentioned in the Introduction, it is natural to expect condensation in real space to 
saturate at some large but finite value for the mass of the condensate (e.g., in shaken 
granular gases when the finite size of the compartments prevents true condensation). 
We will therefore consider systems where the hop rate per particle v{n) asymptotically 
decreases to a finite but non-zero value. Even if the total hop rate per site, u{n), initially 
decreases with n, it eventually starts increasing again and there is no phase transition, 
as discussed above. However, we will show that the initial decrease in v{n) may suffice 
to destabilize a homogeneous state, whereas the asymptotic growth of u{n) prevents the 
formation of condensates with a mass that diverges with the system size. One thus ends 
up with a steady-state containing an extensive number of finite-sized condensates. Since 
there is no thermodynamic transition, canonical and grand canonical ensembles remain 
equivalent. Therefore, we can perform our analysis solely within the latter ensemble, 
and dispense with the cumbersome constraint on the total number of particles. On the 
other hand, simulations are most straightforwardly conducted in the canonical ensemble. 
Since we shall use both ensembles in the following, we will refer equivalently to free 
energy or grand potential with the understood assumption that chemical potential 
and number of particles N are adjusted so that '^nnp{n\fi) = N. 

For the system to start forming condensates (i.e., high-density sites in a sense to be 
defined more formally below) one needs the single site free energy / (n) to be non-convex, 
that is f"{n) < for some range of n. Flat profiles with such occupancies would then 
be unstable under the dynamics and undergo spinodal decomposition. Treating n as 
a continuous variable, we have f{n) ~ J^\nu{n')dn', and hence the second derivative 
with respect to n is given by 



which implies spinodal decomposition for occupancies such that u'{n) < 0. This is 
exactly equivalent to the condition [17] that the hop rate per particle, v{n) — u{n)/n, 



occur. 



fin) 



u'{n) 




u{n) 
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Figure 2. Left: The grand potential density per site, F{fj,,n), for the choice 
u{n) — 6xp(— A(/)arctan(n/0)) discussed in the text, with vq = 2.5, A = 0.01, 
(f) = 250 and (n) = 100. The region unstable to spinodal decomposition is in blue: 
it corresponds to the concave part of the grand potential. Right: The resulting, 
normalized, probability distribution. 



should satisfy the equation 

// N v(n) , , 

v\n) < . 12 
n 

Such an instabihty could in principle lead to complete condensation: the criterion 
(12) is indeed satisfied when the condensation transition occurs. This is because the 
phase separation can lead to phase coexistence between a low-density phase and a high- 
density phase whose density can diverge with system size, i.e., a macroscopic condensate. 
This is for instance what happens in [13], where the jump rates u{n) depends on the 
system size L in such a way that the mass of each of the multiple condensates diverges 
in the thermodynamic limit. For the high density sites to have finite occupancies, we 
will further require that f'{n) — )■ oo when n oo. Combined with the fact f{n) is 
not everywhere convex this implies that the free energy per site has a double tangent 
between two finite densities [18]. Under these conditions the grand canonical potential 
per site, F{n, fi), forms a double well whose minima occur for finite values of n and give 
the typical occupancy of the high- and low-density sites. 'Saturated condensation' now 
is to be understood as referring to this scenario. In terms of the microscopic jump rates, 
the requirement that f'{n) — )■ oo is equivalent to u{n) — )■ oo as n — )• cxo, so u{n) must be 
an unbounded increasing function of n or, correspondingly, v{n) must either increase as 
n — )■ oo, or decrease more slowly than 1/n. 

Let us illustrate with a concrete example of the function u{n) that leads to saturated 
condensation, and which we will use repeatedly throughout the remainder of this work. 
It reads 

u{n) = VqU exp[— A0arctan(n/0)] . (13) 

With this definition, the hop rate per particle v{n) initially decreases exponentially with 
n before it saturates at a constant value 6xp(— A07r/2). A fiat profile is unstable if 
u'{n) < 0, that is if > 2 and n G [p~, p"^] where 

P = 7^ • (14) 
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As n — 7- oo, u{n) ~ fonexp(— A</)7r/2) and is thus increasing linearly with n\ the 
condensates have finite size. This choice of rates gives rise to the double-well free 
energy and bimodal probability distribution shown in fig. 2. 

Note that the range corresponds to the concave part of the free energy, 

as expected from standard thermodynamics [19]. As long as the average density lies 
between the two minima of n), the steady state will be dominated by configurations 
with condensates. Flat profiles will however be metastable outside [p~,p+], thus 
requiring activated events to lead to condensation. Simulations of the systems for rates 
obeying (13) show the predicted behaviour. On the left panel of figure 3 one sees the 
results of simulations started with particles distributed randomly over the L sites. 
The average density p = N/ L is chosen either within or outside the condensation regime. 
The criteria for the condensate to have finite mass can be checked by considering the 
family of rates defined by 

Va{n) = t>on~"exp[— A0arctan(?2/0)] (15) 

One indeed sees that the minimum in the grand potential corresponding to the high- 
density phase is at a finite value of n for a < 1 and diverges when a — )■ 1 (see right 
panel of figure 3). 




Figure 3. Left: Steady-states of stochastic simulations for cf) — 250, A = 0.01, 
Wo — 2.5. For these parameters, the minima in the grand potential correspond to 
PL ~ 60 and pH ~ 800. Occupancies are averaged over a time window t E [5000; 15000]. 
Green and blue symbols correspond to initial average densities (p = 20 and p — 900 
respectively) that are either below p^ or above pH and are as expected stable. For an 
initial density p — 300 between p^ and pn (red symbols), steady-state configurations 
typically exhibits a low density background at p = pL and high density condensates 
at p = ph- Right: Semi-log plot of the typical mass of the high density phase for 
different value of a, using the rates (15). One sees that when a -H- 1, the mass of the 
condensate diverges as expected. 

2.3. Condensation and evaporation dynamics in the steady state 

When the average density lies between the two minima of the grand potential, there is 
a coexistence of high- and low-density sites in steady-state. Since the grand potential 
barrier between them is finite (see figure 2), the instantaneous number of condensates 
will fiuctuate as low-density sites condense and high-density sites evaporate. To discuss 
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these processes in more detail, it is helpful to formally define a condensate (or high- 
density site) as a site with a density greater than the one at the peak in the grand 
potential between the two minima. Likewise, when the density on a site lies below the 
value corresponding to this peak, we refer it as a low-density site. With this definition, 
the average densities of the high- and low-density sites in the steady state are 

while the average number of condensates in steady state is given by 



^-^-L 5] p(n|/x) = L^^. (17) 

_ PH - PL 

— '^peak 

We may then also define the rate of evaporation of the condensates, T^-evap, and 
condensation of the low-density sites, 7?.cond- At equilibrium these balance in such 
a way that the number and size of the condensates remain constant on average: 

7^evapr^^'^=(i^-r^^^)7^cond. 

Although formation and evaporation of a condensate requires many hops and are 
thus complicated processes, the corresponding rates can be computed following a first- 
passage time approach [20]. Indeed if we know the first passage time Tn^np^ak froi^ ^ 
high density site n pn to n — ripeak, we can approximate the evaporation rate for a 
site with n particles by 

7^evap(r^) = ^ , (18) 

where the factor of a half arises as a site at the peak can fall in either direction with 
equal probability. To calculate the evaporation and condensation rates in the steady 
state, we must therefore calculate the relevant first-passage times for diffusion in a 
double well. To achieve this wc follow previous approaches applied to ZRPs undergoing 
a thermodynamic condensation transition [12, 21]. 

We illustrate this procedure by computing the first-passage time to an evaporation 
event. We thus consider a high-density site with n > npeak particles. The rates at which 
the occupancy decreases or increases are given by 

W{n ^n-l)= u{n) = 2 n e-^'^'^^'^M?) (19) 

W{n -^n+1) =ul . (20) 

In principle, the rate at which particles are added on top of a condensate depends on 
the neighbouring densities. Wc will however assume their fiuctuations to be small and 
consider -u^^ to be constant. The first-passage time from n particles to npeak is denoted 
^n,npeak smd, in continuous time, is the solution to the equation 

7;,npeau = + [1 - {u^ + u{n)) dt ] r„,„^^^, (21) 

+ U{n) dt T;-l,np,,k + '^L dt Tn+l.npeak 
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Figure 4. Semi-log plot of the cumulative distribution function of evaporation or 
condensation events in steady state for (f) = 50, A = 2.8/50, wo — 2.5, (n) = 100. The 
blue line corresponds to the cumulative distribution function of a Poisson distribution, 
FcD{t) = exp(— 7t) with 7 « 0.0036 as predicted by equation 24. The red dots stem 
from 10000 simulations and can be fitted with a rate 7 « 0.0031 (green line). 



Lengthy but standard algebra (see Appendix A) leads to 

n ^00 

^-Pe.= E ^I(7)^E^(-^I/^)- (22) 

'="pcak + l ^ '^^ m=l 

Similarly one finds for the first passage time for condensation 

^peak ^ 7—1 

r;w= E (23) 

i=n+l " j=0 

Note that to determine more accurately the rates of evaporation and condensation 
we should in principle average over all starting positions above and below the barrier, 
respectively: 

_ ^n=«pcak P(^\^) ^"."pcak _ Y.n=t Pj^lf^) ^nnp,,k /r, 

We can then define the escape rate from a configuration as the rate for either an 
evaporation or condensation to occur: 

Tie L — Uc , . 

/vtotal = WTf, h —7= . (25) 

evap ^ ^ cond 

In this picture, the distribution of times between events, either evaporations or 
condensations, for a system of length L will be Poissonian with rate T^totai- For the 
choice of parameters = 50, A = 2.8, vq = 2.5, L = 5000 and (n) = 100, equation (25) 
can be evaluated numerically. First we compute n so that (n) = 100 by solving (10). We 
can then use expression (9) for p(n|/i) to compute Tg^ap and Tcond from (22), (23) and 
(24). We last obtain from (25) that the total rate is T^-totai ~ 0.0036. To compare this 
theoretical prediction with numerics we compute the cumulative distribution function 
of evaporation and condensation events from 10000 runs: Fcoit) is the probability 
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Figure 5. Average density in the high density sites as a function of time for parameters 
(f) — 250, A — 0.01, vq — 2.5. The early dynamics {t ^ 0—1000) see the rapid formation 
of Uc condensates that are rapidly growing. The late stage dynamics t ^ 10'^ — 10^ 
correspond to formation and evaporation of condensates that leads Uc to n°!^ and the 
average mass of the condensates to its equilibrium value. Steps in the average density 
correspond to the evaporation of a condensate that is redistributed on the surviving 
ones. (Note the switch from linear to logarithmic scale on the time axis at t — 1000.) 



that the first evaporation or condensation occurs after time t. The simulation data are 
shown in fig. 4; as predicted by the theory the distribution is Poissonian and a fit to the 
simulation data gives 7?.totai ~ 0.0031, which is within 10% of the predicted value. 

3. Two-stage dynamics of condensate formation 

We now examine the relaxation of the system to its steady-state, which is a nontrivial 
process. Starting from a homogeneous configuration within the unstable region {n G 
[p^,p^], see figure 2), the dynamics divides naturally into two regimes presented on 
figure 5. The early-time dynamics see the instability of the fiat profile give birth to 
some number Uc of condensates that then rapidly grow. At the end of this growth 
stage, Uc is in general not equal to the equilibrium number n^^^ and the system has 
not yet reached stationarity. A second stage then follows, taking place on much longer 
timescales, during which activated events responsible for condensation and evaporation 
of condensates lead the system towards its ultimate steady state. This difference in 
relevant timescales between the two stages can be seen from figure 5. Our aim in this 
section is to understand these two distinct relaxational regimes. 

Simulations of the stochastic system, started from random deposition of N particles 
over the L sites of the lattice with N/ L G [p~ , p"*"] , show that the positions of the 
condensates are, initially, anti-correlated (see figure 6). Transient anti-correlations of 
this type are a general feature of systems obeying a conservation law, see, for example, 
reference [22]. This can be understood as a consequence of the condensates being 
created through depletion of the neighbouring sites, thereby preventing the formation 
of other condensates in their immediate surroundings: if there is a condensate at site i, 
there is a decrease in the probability to find another condensate in its vicinity. These 
correlations survive until the late-stage dynamics when new condensates are formed and 
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Figure 6. Snapshots of the correlation function gi{t) = {nj+i{t)nj{t)) ^ / {^nj{t)') ^, 
where {x'^)c = {x'^) ~ (x)'^ and the averages are taken both over the lattice site 
J and many simulations. Starting from an initially flat profile, an anti-correlation 
between sites forms as the condensates condense which then gradually disappears at 
late times when subsequent evaporations and condensations randomize the positions 
of the condensates. The parameters of the simulation are vq = 2.5, $ = 50, A = .05 
and (n) = 80. 



old condensates evaporate, thus smoothing out the correlations |. Note that starting 
with the correct number of regularly- spaced condensates at the correct steady-state 
density would lead to smoothing of the correlations on the same timescales: correlations 
are mainly due to the immobility of the condensates, which only get randomized (by 
evaporation/condensation) in the late stage of the dynamics. Also, if one starts with 
a global density within [pl, Ph] but well outside [p~ , p'^], there is no initial instability 
since the flat profile is metastable and the first stage is thus absent: the creation of 
condensates is then only due to activated events. We now turn to a more detailed 
analysis of both stages. 

3.1. Initial instability and growth stage 

The early stage dynamics corresponds to the rapid growth of an instability around the 
fiat profile which leads to the formation of some number of condensates ric. These 
condensates then rapidly grow and saturate at a density that generically differs from 
that at the minimum of the high density well in the free energy. Insight into this part 
of the dynamics can be gained by comparing the stochastic dynamics of the system 
with its deterministic mean-field limit. The latter is obtained by replacing {u{ni)) by 
u{{ni)) = u{pi); it reads 

Pi = u{pi^i) + - 2 u{pi) . (26) 

and can be integrated numerically using, e.g. , a simple Euler scheme. Starting from an 
initial condition obtained by distributing at random particles among the L sites of 
the lattice, we see in figure 7 that stochastic and mean-field dynamics agree very well, 
despite the fact that the mean-field approximation (by definition) neglects both noise 
and correlations. We infer from this that activated events and spatial correlations are 
not very important to understand the early-stage dynamics and we shall thus proceed 
using this more analytically-tractable mean-field framework. 

X because the model is factorable there can be no correlations in the true steady-state 
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Figure 7. The evolution of the average high and low densities: using numerical 
simulations in the deterministic (red) and stochastic (blue) cases. Both simulations 
were run with random (Poissonian) initial conditions and with the parameters 4> = 250, 
A = 0.01, Wo = 2.5 and (n) ~ 130. Although the agreement is not exact the qualitative 
behaviour is certainly similar. The slight lag between the stochastic and deterministic 
cases is due to activated events increasing the initial separation between high and low 
density sites and is not especially relevant to an understanding of the dynamics. 




At early times, the number of condensates that are created depends strongly on 
the initial condition. For instance, starting the mean-field simulation from a fiat profile 
superposed by a cosine wave leads to the creation of a condensate from each peak of 
the cosine wave, as shown in the left panel of figure 8. We shall first focus on the case 
where the number of condensates is thus controlled. As time goes on, the mass for the 
condensates is drawn from neighbouring sites, which suggests a model of this particular 
condensate-formation process as one in which sites have one of two time-dependent 
densities. Specifically each high-density site has a density pnit), and is surrounded by 
a pair of low-density sites, both with density Phif)- 

Within this simplification, the mean-field equations (26) reduce to 

= 2 [u{pl) - u{ph)] ; Pl = Pec [uipn) - u{pl)] (27) 

where we used the fact that condensates have two low density neighbours whereas low 
density sites have a probability 

Pic = (28) 

of being next to a condensate. This simple approximation reproduces surprisingly well 
the mean-field dynamics for this particular initial condition, as can be seen from the 
right panel of figure 8. Inspection of equation (27) then gives a simple picture of what 
is happening: in the spinodal region u'{n) < 0, so that u{pl) — u{ph) > and pn is 
positive whereas pi is negative. Consequently the high density will increase and the 
low density decrease. This continues as long as u{ph) < u{pl) but stops at the first 
moment when u{ph) = u{pl). This is indeed what happens during the simulation, as 
can be seen on figure 9. 
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Figure 8. Left: Starting from an initial profile which is flat plus a cosine perturbation 
condensates grow at every peak [ric is 35% higher than n°^). Increasing time is 
represented by a change in colour. Parameters are A = O.Of, = 250, vq = 2.5 
and (n) = f30. Right: The mean-field stochastic simulations (red), as left, and the 
solutions to equations (27) (blue). The approximation to consider just two densities 
gives a sharper change but the end points are in excellent agreement. 




Figure 9. Snapshots of the hopping rate out of a site with n particles, u{n) = 
Wo nexp(— A(?!)arctan(n/(}!))), for mean-field simulations (crosses) and stochastic ones 
(circles) . The red and blue symbols represents the average high density and low density 
sites. The blue lines show the steady state values of pH and and the black line the 
function u{n). At t — 400, u{ph) — u{pl) and the condensate thus stop increasing. 
One must then wait for activated events in the stochastic simulations to get closer to 
the equilibrium values. The difference between stochastic and mean-field comes from 
the different number of condensates that results from the initial instability. 

Note that according to the previous discussion, the average values of pi and pn 
at the end of the growth stage can be deduced from the initial number of condensates, 
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Figure 10. Typical evolution of the average density and number of condensate sites. 
Each steps correspond to the evaporation of a condensate that is redistributed over 
the other high density sites. The parameters of the simulations are A = 0.01, = 250 
and vq — 2.5 



using conservation of mass and requiring that u{ph) = u{pl) (see equation (27)) 

UcPh + (L - nc)pL = N ; u{ph) = u{pl) ■ (29) 

Starting simulations with different wavelength for the initial cosine perturbation indeed 
leads to density pn and pl predicted by (29) where Uc equals the number of peaks of 
the cosine wave. 

For more general initial conditions or in the stochastic case, we have not been able 
to find a simple way to predict the number of condensates to be formed. The saturation 
of their growth once formed however follows the same rules as above, leading to high 
and low densities that in general differ from the steady-state ones (see figure 9). Further 
changes in the average densities are due to activated events which change the number 
of condensates and increase both pi and ph- These events are not captured by the 
mean-field approximation, and require a distinct analysis that is discussed in the next 
section. 

3.2. Activated events and late-stage dynamics 

To understand the late stage dynamics, which is mediated by stochastic nucleation and 
evaporation of condensates, we must investigate the activated crossing between the two 
wells in the grand-potential landscape shown in figure 2. 

We have observed that during the late-stage relaxation, the total mass in the high- 
density sites remains approximately constant, despite the fact that the number of such 
sites changes over time. Thus, the average density per condensate increases in the step- 
wise fashion depicted in figure 10. When a condensate evaporates the excess density 
on that site is redistributed over the other high-density sites, thereby increasing their 
average density; conversely, when a new condensate forms, the average density decreases. 
Depending on whether the late stage starts with too many or too few condensates, 
evaporation or condensation will first dominate, before the two rates become closer and 
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Figure 11. Left: the free energy per site for VQ—2.b, (f) — 50, A ~ 3/50, (n) = 100 
at steady state, where there are 40% of condensate sites. Right: the free energy for 
the same microscopic rates and average occupancy but constrained to have 60% of 
condensate sites. 



closer. Once evaporation and condensation of new condensates balance, the steady state 
described in section 2.3 is reached. 

In the steady state, the rates of evaporation and condensation can be treated as 
a first passage problem in a grand potential landscape — as shown previously in section 
2.3. Away from steady state, however, the rates are different from the equilibrium ones 
and depend in general on the dynamics. Numerically, we can measure the rates as a 
function of the fraction of condensates by recording how long the system spends in a 
given configuration with condensates and averaging over many runs. We now show 
that these nonequilihrium rates can also be calculated by appealing to a fluctuation- 
dissipation-type argument and adapting the equilibrium formalism correspondingly. The 
key to this approach is to assume that when the number of condensates ric is sufficiently 
close to its equilibrium value rf^, the difference — rf^ could be due either to a 
spontaneous fluctuation (which we are observing) or equivalently to the application of 
a small field. To this end, we introduce a new 'doubly-grand canonical' ensemble that 
involves an additional chemical potential yu' conjugate to the number of condensates ric- 
The corresponding partition function is then given by 

oo L 

2 = 5^ 5^e'^^e^'"^Zi,7v . (30) 

Ar=0 nc=0 

Defining a condensate, as before, as any site containing more than npeak particles, the 
partition function can be re-written as 



.n=0 



-f{n)+^n 



E 



E 



exp[ 



-G(/i, /i', n)] 



.n=0 



(31) 



where we have introduce a new thermodynamic potential 

G(/i, /i', n) = f{n) - nn- n' 9{n- ripcak) (32) 

G{fi,fi',n) can simply be obtained by introducing a step — /i' in the grand potential 
F{fi,n), see fig. 11. The new marginal probability to observe an occupancy n is then 
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Figure 12. Red points show data from stochastic simulations of the microscopic 
dynamics. The time spent in a configuration with ric condensates was averaged over 
multiple runs and the total rate to leave a given configuration was taken as the inverse 
of this average time. The blue curve shows the rate calculated by using equations (24) 
to determine the first-passage time to the peak and assuming a site at the peak is 
equally likely to fall into either well. All data calculated with vq = 2.5, (p = 50, 
A — 3/50, {n) — 100 and on a lattice of length L = 500 sites. 



given by 

— G(p n' n) °° 

p{n\^^,^^) = %——■, where Hi(/x, /i') = e-^^'^''^' (33) 

To evaluate the rates of evaporation and condensation of the system in the presence of 
He condensates, one can thus compute numerically the values of /x and /i' such that 

^np(n|/i,/x') = and 9{n - npeak) p(^|/i, jj!) = ng (34) 

n n 

We can then compute the rate of evaporation or condensation from a configuration with 
N particles and Uc condensates exactly as in section 2.3 where F{fi,n) is now replaced 
by G(/i, /i', n). To compare the results of these calculations with numerics we started 100 
runs from Uc = 275 and 100 from Uc = 155. While these runs relaxed to the equilibrium 
value ric = 199 we recorded the average time spent, r(nc), by the system for each 
intermediate value of Uc and approximated 7?.totai(^c) = 1/T(nc). The results of these 
simulations are compared with the theoretical predictions in figure 12. The fact that the 
agreement between theory and numerics is not as good as in section 2.3 may be due to 
the fact that we have poorer statistics (100 runs against 10000). Nevertheless, we find 
that this 'doubly- grand canonical' construction provides remarkably good estimates for 
the evaporation and condensation rates: even where they are two orders of magnitude 
larger than in equilibrium, theory and simulations are still within a factor 2 of each 
other. 
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4. Conclusion 

In this work we have identified a saturated condensation scenario that may occur within 
the zero-range process (ZRP). The steady state is characterized by an extensive number 
of finite-sized condensates, as opposed to a single macroscopic condensate that has 
previously been the focus of attention. Such a state may be brought about, for example, 
by a constraint on the total mass that may occupy a single site. We have determined 
the conditions on the hop rates within the ZRP that must be satisfied for saturated 
condensation to arise, and have investigated various aspects of the model's dynamics, 
both within and en-routc to the steady state. 

Our results complement existing treatments of the condensation dynamics for 
the version of the ZRP with a true condensation transition [9, 10, 11, 12, 13]. In 
those works the focus has mostly been on the late-time scaling behaviour of cluster 
size, and the characteristic timescales for evaporation and reformation of condensates. 
Within the model of saturated condensation discussed here, we have obtained a more 
complete account of the dynamics — including the nontrivial behaviour of the number of 
condensate sites — from very early times right through to the steady state. In particular, 
we find that the relaxation takes place in two stages: first, some number of condensate 
sites is dynamically selected which depends on the initial condition. These condense 
rapidly, but leave the system in an out-of-equilibrium state that slowly relaxes through 
activated evaporation and condensation events. A mean-field approximation proved 
reliable in analyzing the first stage, and a first-passage calculation conducted within a 
specially-constructed grand canonical ensemble well described the second. 

There are, however, some discrepancies between these approximate analytical 
approaches and the full stochastic dynamics. For example, the rates we calculated for 
evaporation and condensation, both at equilibrium and away from the steady state, seem 
to overestimate those measured from the simulations. One plausible explanation for this 
inconsistency is the assumption that the neighbouring sites will have a constant density 
throughout the evaporation or condensation process. In fact the neighbouring sites are 
likely to have higher densities than otherwise during an evaporation and lower densities 
during a condensation. As both high and low density sites sit in regions where v'{n) > 0, 
an increase in density will increase the rate at which particles enter the evaporating site 
whilst a decrease in neighbouring density will decrease the rate at which particles enter 
a condensing site. Both these processes will have the effect of increasing the average 
duration of a condensation or evaporation event, and reducing the rates at which they 
happen. Hence the calculated evaporation and condensation rates would be larger than 
the real rates as measured. 

More serious, perhaps, is the absence of a satisfactory explanation for the number 
of condensate sites formed from a uniform initial condition. One would anticipate that 
the early-time noise would play a major role in determining this number, although we 
have been unable to relate these two quantities directly. This we leave as a possible 
avenue for future work. 
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Finally, since we believe our model to reproduce more faithfully experimental 
situations such as the shaken granular gases than the traditional version of the ZRP, 
it would be very interesting to investigate both early- and late-time dynamics of the 
corresponding experiments [2]. 
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Appendix A. Solution of the first paissage time problem 

The first passage time from n particles to ripeak is denoted T^^npe^i^ and, in continuous 
time, is given by tfie solution to the equation 

^n,np,,k = dt + [1 - (ml + U{n)) dt ] Tn-np,^^ + u{n) dt Tn-l,np,,^ + Ul dt T„+l,„p,,k (A-1) 

This equation states that the time to go from n to npcak particles (l.h.s) is dt plus the 
time to go from the new number of particles, obtained after a time interval rft, to npeak- 
With probabilities u{n) dt and ul dt. there arc now n — 1 or n + 1 particles, while with 
probability 1 — {ul + u{n)) dt there arc still n particles, hence the three terms of the 
r.h.s. Equation (A-1) then reduces to 

with boundary condition Tn^^^^^npea.k — 0- Now, define the difference dn = Tn^^peak ~ 
^n-i,npeak so that 

u{n) dn - Ul dn+i = 1. (A-3) 
Note that the evolution of the probability to find n particles at a site, p{n\ii), is given 

by 

— — — = u{n + l)p{n + l|/x) + ULp{n - - {u{n) + ml) p(n|/x) (A-4) 

= Jn+l,n Jn,n—l (A"5) 

where Jn+i,n = ^{n + l)p{n + l|;u) — ULp{n\^). At steady state the left hand side 
of this equation must equal zero, so the current, J must be constant. However, as 
p{n < 0|/i) = and Vq = 0, J-i^ = and hence this constant must be zero. This 
implies that u{n-\- l)p{n-\- — ULp{n\ii) so the solution to the homogeneous version 
of equation (A-3) is: 

dn = ( \ \ \ - (A-6) 
u[n)p[n\ii) 

To solve the inhomogeneous equation, then, we look for solutions of the form 

dn^ , I (A-7) 

u\n) p{n\jj,) 

where c„ is to be determined. Substituting this expression back into equation (A-3) 
gives 

—Ul —, + uiri) , , , = 1. (A-8) 

^^/(n + l)p(n + l|/x) ^ ^ u{n)p{n\ii) ^ ' 

The detailed balance condition, u{n -\- l)p(n + l|/i) = ULp(n|/x), then implies 

Cn - c„+i = p(n|/i), (A-9) 
so that Cn = Si^„p(^|A*) and the first passage time is given by 

= E syy^EM^lM). (A-io) 



T 



